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Abstract 

Hypersonic laminar flow over a quilted surface 
configuration was Investigated to determine the 
pressure and heating loads to the surface. The 
quilted surface configuration simulates an array of 
metallic Thermal Protection System panels that are 
bowed In a spherical shape as a result of thermal 
gradients through the panel thickness. Panel 
bowed heights" (deflection of panel from its 
original flat configuration) from zero to 1.25 
times the local boundary layer thickness were 
Investigated for flows both longitudinal ly and 
diagonally aligned with the panels. The flow field 
over the configuration was mathematically modeled 
using the time-dependant three-dimensional 
conservation of mass, momentum (Navler-Stokes), and 
energy equations. A boundary mapping technique was 
used to obtain a rectangular parallelpiped 
computational domain, and a MacCormack explicit 
time-split predictor-corrector finite difference 
algorithm was used to obtain steady-state 
solutions. Results show that for bowed heights 
less than or equal to the local boundary layer 
thickness total Integrated heating loads vary 
linearly with bowed height. The total Integrated 
heating load on the leading panel and the 
downstream panels Increase to less than 1.44 and 
1.26 times the flat plate value, respectively. 

When the bowed height exceeds the flat plate 
boundary layer thickness heating rates Increase 
more rapidly, vorticity effects become significant 
and flow begins to separate ahead of the leading 
panel. 
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Subscripts 

max 

cn 


functional relationship between 
the physical and computational 
domains (eqs. A7 and A9) 
coordinate directions 
flat plate boundary-layer 
thickness at x =» xm 
one-half the fillet width 
transformed coordinates 
redistributed transformed 
coordinate (eq. A8) 
viscosity (eq. 6) 
density 

stress tensor (eq. 2) 


maximum coordinate value 
coordinate value at center of 
panel n (n « 1, 2, 3) 


Introduction 


Metallic Thermal Protection Systems (TPS) have 
been proposed as viable concepts for hypersonic 
vehicles and as an alternative for the current 
Reusable Surface Insulation (RSI) TPS on the Space 
Shuttle. Metallic TPS would be less fragile 
than RSI, last longer, eliminate gaps, and offer 
the advantage of mechanical fastening. Concepts 
combining thin foils of titanium and/or superalloy 
materials sometimes combined with metallic 
honeycomb and/or Insulation materials have been 
fabricated. These concepts are weight competitive 
with RSI. 
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Nomenclature 


specific heat at constant volume 
(eq. 5) 

internal energy 

total energy (e + 1/2 (u^+v^+w^)) 
bowed panel height above flat 
plate 

thermal conductivity of air 
variable grid concentration 
parameter (eq. A8) 
differential operators (eq. Bl) 
pressure 

gas constant (eq. 4) 

radius measured from panel center 

radius of flat-plate/spherical 

region Intersection 

temperature 

time 

velocity components In x, y, and z 
directions, respectively 


However, the use of metallic TPS can Introduce 
some physical phenomena not present In the RSI 
system. One such phenomenon Is thermal bowing of 
the metallic TPS panels.^ Bowing occurs because 
metals have relatively large coefficients of 
thermal expansion and because In the course of 
performing their function as a thermal protection 
system the TPS panels must sustain large 
temperature gradients through their thickness. 

These conditions lead to larger thermal expansion 
of the outer portion of the panel than of the Inner 
portion. The square panels are held down at the 
corners but otherwise thermal expansion Is 
unrestrained to minimize thermal stresses. As a 
result the center of the panel bows out Into the 
stream so that the vehicle moldline Is altered. An 
array of panels would take on a lightly quilted 
configuration. 

High speed flow over a thermally bowed 
metallic TPS Introduces complicated 
three-dimensional flow phenomena that perturb the 
flow field and significantly alter the wall 
pressure and heating rate distributions. No 
analytical studies of a configuration resembling 
the thermally bowed TPS were found In the 
literature. A preliminary study of flow over 
spherical dome protuberences showed that the domes 







did significantly effect the flow fields. ^ in some 
cases separated flow regions developed and trailing 
vortices formed, significantly altering downstream 
heating rates. 

This study was conducted to determine pressure 
and heating rate distributions on a quilted surface 
configuration simulating an array of thermally 
bowed metallic TPS panels In a high speed laminar 
flow field. The flow field simulates Mach 7 flight 
at 120,000 ft. altitude and is similar to the test 
condition In reference 4. Sections of arrays two 
and one-half panels In length have been modeled for 
flows longitudinally and diagonally aligned with 
the panels. Computer storage limitations 
restricted the model to 0.31 ft. square panels. 
Bowed panel heights, h, in the range 0 £ h £ 0.25 
inches were considered In a flow field whose flat 
plate boundary layer height, 4fp, Is 0.20 Inches. 
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Ideal gas state equation. 


The analysis employs the same formulation and 
numerical procedures used In the preliminary study 
of spherical dome protuberances^ a three-dimensional 
formulation of the time dependent equations for 
conservation of mass, momentum (Navler-Stokes), and 
energy and a MacCormack explicit time-split 
predictor-corrector finite difference algorithm. 
Steady-state solutions were obtained on Langley's 
Cyber 203 vector processing computer. The results 
of those solutions are presented In this paper using 
state-of-the-art computer graphics techniques. 


P “ pRgT (4) 

These equations along with the Ideal gas 
relationship, 

e = CyT ( 5 ) 

and the Sutherland viscosity formula, 

3/2 

“ “ 2.270 X 10-8 (6) 


Analysis 


Governing Equations 


Air, the fluid medium, was modeled as a 
compressible, viscous, thermally conducting gas and 
was assumed to be an Isotropic newtonlan Ideal gas 
In local thermodynamic equilibrium. The 
three-dimensional time-dependent system of 
equations describing this model Include the 
continuity equation (conservation of mass), the 
Navler-Stokes equations of motion (conservation of 
momentum), the energy equation (conservation of 
energy), and the Ideal gas equation (thermodynamic 
state equation). These equations have been derived 
In the literature” and are shown below written In 
cartesian tensor form where repeated Indices Indi- 
cate a summation over the range of Indices (I.e. 
3(PWJ ; 
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3X 


aXj ^ 3X2 3X3 


where Wi, Wg, Wg are the velocity components u, v, 
w, and A 2 » A 3 are the coordinate directions x, 
y, z. The Kronecker delta* equals 1 when 1 » 
j and equals zero when 1 / j. 


Continuity equation^ 
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Navler-Stokes equations of motion^ 
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constitute a set of five partial differential 
equations and one algebraic equation In six 
dependant variables: three velocity conponents, u, 
V, and w; density, p; temperature, T; and pressure, 

p • 


The physical domain grids and array dimensions 
for the quilted surface configurations with 
longitudinally and diagonally aligned flow are 
shown In figure 1 . Points corresponding to corners 
of the metallic TPS panels were fixed while the 
remainder of the panels bowed In a spherical manner 
simulating the thermal expansion that would occur 
during aerodynamic heating. The resulting surface 
topography Is a lightly quilted pattern. Leading 
panels were faired Into the upstream flat plate by 
extending the spherical shape. The Intersections 
between the spherical regions and the flat plates 
were made mathematically smooth by modeling a 
fillet region at '*0 ± e with a Hermite polynomial 
(where rQ = 0.21 ft. and e » 0.01 ft.). Panel 
"bowed heights" (deflection of panel from original 
flat configuration) between 0 £ h/«f„ < 1.25 were 
modeled In 0.25 h/fi^p Increments. The 
longitudinally aligned flow conHguratlon has a 
longitudinal symmetry plane at the center of the 
panels and another at the edge of the panels, 
reducing the required computational model to 
one-half the panel width (fig. la). The diagonally 
aligned flow configuration has longitudinal 
symmetry planes on the panel diagonals, reducing 
the required computational model to one-half a 
panel diagonal (fig. lb). The domains of the 
longitudinally and diagonally aligned flow 
configurations were discretized Into 94 x 32 x 26 
and 76 X 45 X 26 grids, respectively. The 
resulting computations contained approximately 5 x 
10 ® degrees of freedom. 
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Boundary Conditions 

A boundary layer profile was Input as the 
upstream boundary plane (x ■ 0) data and held 
constant throughout the confutations. The boundary 
layer profile was determined using a two-dimensional 
boundary layer analysis' for flow over a flat 
plate. The length of flat plate ahead of the three 
dimensional region was selected to produce a flat 
plate boundary layer thickness, = 0.2 In., at 
the center of the leading panel (X^» 0.31 ft.). The 
free stream conditions for Mach 7 flight at an 
altitude of 120,000 ft. used to Initiate the 
calculations are shown In Table 1. Using the 
converged flow field solution for one bowed height 
as the Initial condition for the next greater bowed 
height reduced the time required for convergence bv 
approximately one-third. ^ ^ 

K 4 boundaries are planes of symmetry on 

which the dependent variables are determined by 
quadratic extrapolation from Interior points so 
that the cross flow gradients, 3/3z, of all state 
variables are zero and the cross flow velocity, w. 

Is zero* The upper boundary and downstream 
boundaries are specified as no-change boundaries, 
i.e. , set equal to the plane next to them. The 
wall Is a no-slip boundary with a constant 
temperature of 1440®R. This wall temperature Is 
both the upper use temperature and the flat plate 
radiation equilibrium temperature for the 
Input boundary conditions for a titanium TPS 
assuming a high einittance coating and zero thermal 
diffusivlty. 

Computational Domain 

In order to generate a boundary fitted 
coordinate system and provide Increased resolution 
In the boundary layer, the physical domain (fig. 1) 
was transformed Into a rectangular parallelepiped 
computational domain using a two-boundary grid 
pneratlon technique.^ The computational domain 
has equally spaced grid lines In the C, n, and c 
directions* Transformation of the governing 
equations to the computational domain Is shown In 
appendix A* 

Numerical Procedure 

A finite-difference solution to the governing 
equations was attained using the second-order 
accurate HacCormack explicit time-split 
predictor-corrector algorithm.^ A fourth-order 
eddy-viscos1ty-l Ike damping term was Included to 
promote stability In the regions of large pressure 
gradients. Both the algorithm and the damping 
term are further discussed In appendix B* ^ ^ 

Programming and Computer Processing 

The numerical procedure has been programmed to 
run on Langley s Cyber 203 vector processing 
machine using 32-b1t (half word) arithmetic, to 
provide Increased speed and storage. Vector 
processing and the large data base required an 
Interleaved storage data management system to 
promote efficient processing* The data management 
system and vector processing are described In 
reference 8. 


Results and Discussion 

The effects of hypersonic flow over a quilted 
surface configuration simulating spherically bowed 
metallic TPS panel arrays with both longitudinally 
and diagonally aligned flow were analyzed* Bowed^ 
heights from zero to 1.25 times the local flat 
plate boundary layer height were Investigated. The 
results of the analysis are shown herein with the 
aid of computer graphics displays and vector 
plots. All similar parameters, pressures, heating 
rates, and velocity vectors, are displayed on a 
common scale to facilitate direct comparisons* 


Flow Field Perturbations 

The effects of bowed height on the surface and 
flow field pressure distributions for 
longitudinally and diagonally aligned flow are 
shown In figures 2 and 3, respectively. Conputer 
graphics displays (fig* 2a and 3a) present the 
pressure distributions on the wall, on the two flow 
field symmetry planes, and on the outflow 
plane— over 8400 data points In each display. 
Pressure at the Initial plane Is 0.1614 psi and 
pressure variations range from 0.088 psi to 0.395 
psi for the h/6^ « 1.25 case. Plots of the peak 
pressure for each panel as a function of bowed 
height are shown In figures 2b and 3b* The leading 
panels In both flow configurations generate shock 
waves and have pressure distributions similar to 
those found on spherical dome protuberances ^ 

The peak pressure on the windward side of the 
leading panel Increases exponentially with 
Increasing bowed height to 2.22 times the flat 
plate value at h/6^p = 1.25. 

Downstream panels In longitudinally aligned 
flow (fig. 2) are shielded by the leading panel and 
generate much weaker shock waves than the leading 
panel. The peak pressures on panels two and three 
ncrease linearly with bowed height to 1.33 and 
1.29 times the flat plate value, respectively, at 
"'®fp ■ 1.25 {fig. 2b). Expansion regions form 
on the ]e«1de of each panel resulting In a flow 
field with alternating high and low pressures which 
tend to merge away from the surface. 

Panel two In diagonally aligned flow (fig. 3) 

IS only partially shielded by the leading panel. 
Impinging flow generates a shock wave similar to 
the leading panel but not as strong. The peak 
pressure on panel two Is 1.78 times the flat plate 

Shielded by the leading panel and generates a much 
weaker shock wave than the leading panel. The peak 
pressure on panel three Increases linearly to 1.40 
times the flat plate value at h/6f. » 1.25. Again, 
the alternating compression and expansion regions 
produce a flow field with alternating high and low 
pressure regions. In addition the offset panels 
cause these regions to be Interlaced In the 
transverse direction. 

The effect of surface configuration on the 
flow near the wall for h/«fp = 1.0 Is shown by 
the streamwise/crossflow velocity vectors In figure 
4. Vectors for the longitudinally aligned flow 
(fig. 4a) show moderate three-dimensional flow 
relief on the windward side of the leading dome but 

n«aie the lees Ide or on downstream 

panels. Separated flow regions have developed 
ahead of the lead panel and In the Interpanel 
valleys normal to the flow. Vectors for diagonally 
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aligned flow (fig* 4b) show strong three 
dimensional flow relief on both the windward and 
leeside of all panels. The flow follows the 
diagonal valleys between panels. This Interaction 
of flow field and surface contour generates 
vortices that Impinge on downstream panels. A 
small separated flow region has developed ahead of 
the leading panel similar to the spherical dome 
protuberances in reference 5. In addition, larger 
separated flow regions have developed In the 
depressed corner regions of the panels. Much 
larger separated flow regions develop on the 
windward side of the leading panels for both flow 
alignments when h/«S^p Is increased to 1-25. 

Wall Heating Rate Distributions 

The effects of bowed height on the surface 
heating rates for longitudinally and diagonally 
aligned flow are shown In figures 5 and 6, 
respectively. Computer graphics displays (fig. 5a 
and 6a) present the heating rate distributions on 
the walls. The average heating rate on the flat 
plates is approximately 1.5 Btu/ft^ sec and ranges 
between 0.12 and 4.88 Btu/ft^-sec for h/5^p = 

1.25. Plots of the peak heating rate for each 
panel as a function of bowed height are shown in 
figures 5b and 6b and plots of the total integrated 
panel heat loads as a function of bowed height are 
shown in figures 5c and 6c. Heating rate 
distributions on the windward side of the leading 
panel in both flow configurations are similar to 
those found on spherical dome protuberances, ^ The 
peak heating rates on the leading panels Increase 
exponentially with bowed height to 3.31 times the 
flat plate value at h/^^^p = 1.25. 

Heating rate distributions for longitudinally 
aligned flow (fig. 5a) show increased windward side 
heating on each panel and decreased heating In the 
valleys between panels as bowed height Increases. 

A separated flow region In front of the U/^fp - 
1.25 leading panel slightly perturbs the local 
heating rate. The heating rate distribution on the 
h/<Sf = 1.25 third panel also shows the Incipient 
formation of off-center high heating rate areas 
caused by impinging vortices from upstream panels. 
Peak heating rates on panels two and three Increase 
approximately linearly with bowed height to 2.05 
and 2.02 times the flat plate value, respectively, 
at h/<S^p “ 1.25. Panel three shows a slightly 
non-linear Increase between h/<5fp =1.0 and 1,25 
as a result of increased heating^caused by the 
impinging vortices. Total panel Integrated heating 
loads (fig. 5c) do not Increase as sharply as peak 
heating loads because of reduced leeside heating 
rates. Panel one heat load Increases exponentially 
with bowed height to 1.61 times the flat plate 
value at h/<5fp = 1.25. Panel two heat load 
Increases approximately linearly with bowed height 
to 1.33 times the flat plate value at h/5fp = 1.25. 

Heating rate distributions for diagonally 
aligned flow (fig. 6a) show similar windward side 
and valley heating trends and separated flow 
perturbations. However, panel two and three 
heating rates vary significantly from the 
longitudinally aligned flow configuration as the 
bowed height exceeds h/5^p = 0.50. Panels two 
and three have Increased heating regions on the 
centerline of the panels that are the result of 
flow field impingment- However, Impinging vortices 
also cause regions of higher heating to form 
off-center that result In exponential increases In 


peak panel heating (fig. 6b) to 2.60 and 2.34 times 
the flat plate value, respectively, at h/6^p = 

1.25. Although peak panel heating rates for 
diagonally aligned flow are about the same or 
higher than those for longitudinally aligned flow, 
total panel Integrated heat loads (fig. 6c) are 
less as a result of lower leeside heating. Total 
heat loads on panels one and two Increase 
exponentially to 1.54 and 1.21 times the flat plate 
value, respectively, at h/5^p = 1.25. 

These computations have been done assuming a 
constant wall temperature of 1440®R. However, 
because of the wide range of heating rates 
encountered, the wall temperature will not be 
uniform. The equivalent wall radiation equilibrium 
temperature for the peak heating rate for each 
panel was computed assuming zero thermal 
diffusivity of the wall, a constant heat transfer 
coefficient, and a surface emissivlty of 0.72. The 
resulting peak wall temperatures shown In figure 7 
are conservative. The equivalent wall radiation 
equilibrium temperature for all the oanels 
longitudinally aligned flow (fig. 7a) Increase 
nearly linearly over the range of bowed heights. 
Panels one, tv/o, and three temperatures Increase to 
1826, 1603, and 1573®R, respectively! at h/6fp = 

1.25. The temperatures for the panels diagonally 
aligned flow (fig. 7b) behave In a similar manner 
until the vorticity Induced heating on the 
downstream panels causes a non-linear increase In 
heating. Panels one, two, and three temperatures 
Increase to 1826, 1679, and 1702®R, respectively, 
at h/6^p = 1.25. 

These data show that for bowed panels that are 
submerged In the boundary layer, I.e. , h/6fp < 

1.0, vorticity induced heating is negligible and 
flow does not separate on the windward side of the 
leading panel- Integrated heat loads on the 
leading panels and downstream panels Increase to 
less than 1.44 and 1.26 times the flat plate value, 
respectively. In comparison, under similar flow 
conditions experimental ly measured total heat loads 
on simulated Shuttle RSI tiles Increase an average 
of 1.27 times the flat plate value as a result of 
the gaps between tiles. Local radiation 
equilibrium temperature Increased less than 20 
percent on leading panels and 11 percent on 
downstream panels. 

True spherical deformation of a metallic TPS 
panel, as prescribed In this study, would occur 
only If the heating were uniform. In reality the 
aerothermal heating analysis Is coupled to the TPS 
panel thermal/structural analysis through the 
heating rate distribution. The actual thermal 
deformation will vary from a spherical shape 
because the heating rate Is nonuniform. Solution 
of a fully coupled problem of this magnitude is 
beyond current capabilities. An iterative solution 
between the fluid analysis and the 
thermal /structural analysis would be expensive. 


Conclusions 

An analytical study of high-speed flow over a 
quilted surface configuration simulating 
spherically bowed metallic Thermal Protection 
System panel arrays was conducted to determine the 
associated flow field and the effect on panel 
aerothermal heating. Flow longitudinally and 
diagonally aligned with panels whose bowed heights 
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varied from 0 to 1.25 times the flat plate boundary 
layer thickness was Investigated. The flow field 
was modeled as three-dimensional, laminar viscous, 
and thermally conducting flow utilizing the continu- 
ity equations of mass, momentum, and energy. 

Results show that for bowed heights less than or 
equal to the flat plate boundary layer thickness 
integrated heat loads on leading and downstream 
panels Increase to less than 1.44 and 1.26 times the 
flat plate value, respectively, and radiation equil- 
icnon"’ t®'"P6'"atures Increase by less than 280 and 
loO R, respectively. These Increases are similar to 
those seen on Reusable Surface Insulation tiles due 
to pps. when the bowed height exceeds the flat 
plate boundary layer thickness heating rates 
Increase rapidly, vorticity effects become slgnlfl- 
cant (especially for diagonally aligned flow), and 
flow begins to separate ahead of the leading panel. 
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Appendix A 

Governing Equations In Computational Domain 

Transformation of the governing differential 
equations, eqs. 1 , 2 , and 3 , to the computational 
domain requires replacement of all spatial 

Jo[l^Sg?’ ^ 2 . 3 . Sith the 
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coordinates n snd c, l.e. , 

The resulting 

equations of motion are. 
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Written In short form the system of partial 
differential equations Is, 


3U . 3F 
3t 3C, 

where. 

"p ” 

U = 

PU 

pv 


pw 


-pEJ 


3C. 


3G 


3C. 


+ * X 3H "~t 

3Cjj 3y 3C^ az 


3C, 


(A5) 


5 



F 


pu 



puu 

- m 


puv 

-T12 


puw 

- fi3 


pEu 

V 3T 
- 3C 

m 

3X 


(UTII + VT12 + WTI3) 


where |^, and which vary for the flat 

plate, fillet, and spherical regions of each panel 
n (n « 1, 2, 3) as defined by the radius measured 
from the center of each panel 

'•n " >/ ^'max « ‘ ^ <"max ' ' "cn>^ 
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puv - 1^2 
pvv - t 22 
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PUW - T13 
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where x^n the streamwise coordinate and z-n 
is the spanwise coordinate of each panel center. 

These Jacobian elements are calculated and 
stored for each grid point in the domain thus 
allowing the solutions from the computational 
domain to be transformed to the physical domain. 

Substituting the transform Jacobian (eq. AlO) 
into the short form of the equation system (eq. A5) 
reduces the system to the following, 

at an ax Sn ay %Taz ° 


(A12) 


The required Jacobian transform matrix, 
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is determined by a linear mapping of the boundaries 
at n = 0 and n = 1 of the computation domain onto 
the boundaries of the physical domain at y « 0 and 
y * y^ax the following form; 

Y » Y (c, 1, c) n + Y (5, 0, O (1 - n). (A7) 

Further transformation to increase resolution in 
the boundary layer region entails an exponential 
grid stretching in the n direction of the form. 


n “ 



(A8) 


where k is a free parameter that adjusts grid 
spacing. Equation A7 becomes. 


Appendix B 
Numerical Procedure 


A finite-difference solution to the governing 
equations was attained using the second-order 
accurate MacCormack explicit time-split 
predictor-corrector algorithm.^ The algorithm 
splits the differencing scheme into a series of 
one-dimensional operators. The operator in the 
principle flow direction, Lr, advances the 
computation by one time increment while the 
operators in the other two coordinate directions, 
and L^, are divided into two steps of one-half a 
time Increment each and applied symmetrically 
around the principle flow direction operator. This 
arrangement allows the principle flow direction 
operator to run at a Courant, Friedricks, Lewy 
(CFL) number close to the optimum value of one. 

The operators are applied serially as follows: 


U 


n+1 


‘•c(f) '■5M 



(Bl) 


T “ Y (5. 1, c) n + Y U, 0, c) (1 - n) (A9) 

For the geometry of the subject problem the 
resulting Jacobian matrix Is of the form. 


^max 



In 

in 

in 

3X 

3y 

3Z 

0 

0 

1 

^max 


(AlO) 


Each operator has two steps, a predictor step that 
advances the solution by its time increment based 
on a backward spatial differencing of the 
one-dimensional equations and a corrector step that 
recalculates the advancement based on a forward 
spatial differencing of the predicted results and 
averages the result with the predictor. The 
predictor must therefore lag the corrector by a 
spatial Increment in each coordinate direction. 

The combination of these two steps results in a 
second-order-accurate differencing scheme in time 
and space. As an example the first operation L 
(At/2) applied to the data at time step n is as 
follows: 


6 



Predictor step : 



Corrector step : 



Where bars over the variables Indicate predictor 
values and asterisks over the variables Indicate an 
Intermediate result in the serial application of 
operators as shown In equation Bl. The other 
operators have similar form. 


Although a complete stability analysis of this 
algorithm is not available the CFL limit yields a 
conservative time step as follows. 


At < min 


ImJ + IjlI + IwJ 

Ax Ay Az 


+ c 





(B4) 


where c Is the local speed of sound. 


In the formulation of this algorithm teniis of 
order three and higher have been truncated. These 
terms are not significant In many flow regimes, 
but, when modeling high energy flow fields with 
strong shocks, as In this paper, they can be 
significant. In the region of a shock, where 
pressure gradients are extremely sharp, the 
higher-order terms act as dampers on the pressure 
disturbance. When these dampers are missing, as In 
most finite-difference formulations, they can be 
simulated by adding terms such as the fourth-order 
eddy-viscosity-llke term suggested In reference 10 
which has the form 


a At AC^ 
1 



I u I + c 
4P 



(B5) 


Implementing this smoothing term required adding a 
term of the following form to each vector 
component, both predictor and corrector steps, in 
the F, G, and H vectors defined In equation A5. 


a 


Ei+LM 

’’i+l.j.k 


- 2 P 

TTT 


LM 

1,j,k 


+ P 

TT 




(Vu.k ■ 

where o Is a free adjustable parameter whose range 
is 0 < a < 1/2, c Is the local speed of sound and U 
Is the vector defined In equation A5. The term Is 
of significant magnitude only when the pressure 
gradient Is large. 


Table I Specified conditions 


Free Stream 


Mach Number, M« . . . 
Velocity, u«, ft/sec 
Pressure, P„,lbf/1n2 
Temperature, Teo, R 
Density, p^, lbn,/ft^ 

Wall 

Temperature, T^, R . 
Reynolds Number, Rg . 


7 

. . . 7131 

4.92 X 10 -^ 

.... 432 

4.432 X 10 “^ 


. . . 1440 

3.38 X 1 Q 5 



(a) Longitudinal aligned flow 



(b) Diagonally aligned flow 
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Fig. 1 Physical domain grids and array dimensions 
for quilted surface configurations. 
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(a) Pressure distributions 


PANEL DESIGNATION 


PANEL #1 


PANEL #2 


•PANEL #3 


pressure. 


Pressure on a quilted panel array with 
longitudinally aligned flow as a function 
of bowed height. 
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HEATING RATE, 


ing rate distributions 


PANEL DESIGNATION 


PANEL DESIGNATION 


•PANEL #3 


■PANEL #2 


(b) Peak panel heating rate. 

Fig. 6 Heating on a 
diagonally al 
bowed height. 


(c) Integrated panel heating load 

lel array with 
as a function of 






Ffg, 7 Maximum equivalent wall radiation 

equilibrium temperature as a function of 
bowed heigbt- 
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